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Abstract 

The theory of gravity with an auxihary extra dimension is known to give 
the ghost-free cubic completion of the Fierz-Pauh mass term in the decouphng 
hmit. Our work generahzes the boundary condition in the auxihary dimen- 
sion that avoids ghosts order- by-order, and to ah orders, in the decouphng 
hmit. Furthermore, we extend the formahsm to the case of many auxihary 
dimensions, and we show that the multi-dimensional extension with the ro- 
tationally invariant boundaries of the bulk, is equivalent to the model with a 
single auxiliary dimension. The above constructions require the appropriate 
adjustment of the boundary condition, which we discuss in detail. The other 
possible extension of the original model by the Gauss-Bonnet term is studied 
as well. 



1 Introduction 

The theory of massive gravity which once had been thought to be ruled out because 
of propagation of ghost at non-hnear level [1, 2], was recently revitalized in [3, 4], 
where it was observed that there is a two parameter family of generalizations of Fierz- 
Pauli mass term [5], which is free from ghostlike instability at least in the decoupling 
limit, and also in the full theory at least in the quartic order in nonlinearities [4]. 
For further interesting studies of this model see [6, 7, 8, 9]. This on the other hand 
strongly motivates a quest for underlying setups that would naturally explain or at 
least elegantly reproduce the unusual structures emerged in [3, 4]. 

The theories of gravity supplemented with an auxiliary extra dimension (AED) 
[10, 11], inspired by the DGP brane-world gravity model [12], appeared as a promis- 
ing step in this direction as they correctly predicted a cubic ghost-free completion 
of Fierz-Pauli [13]. In these models while all matter fields live on a 4 dimensional 
brane, the metric has been extended to an extra dimension —1 < -u < 1, and is 
called g^^{x,u). The brane is located at m = and a Z2 symmetry is imposed on 
the fields, g^i,{x,u) = g^u{x, —u), using which the graviton's Lagrangian looks like 



where the first term is the Einstein-Hilbert Lagrangian on the brane as a function 
of g^v{.x) = g^y{x,u = 0), while k^,^ = jdug^u, k = g'^'^k^^, and aU indices in the 
second term are contracted using inverse extended metric g'^^. The coordinate u is 
called an auxiliary dimension because after choosing a second boundary condition, 
say at It = 1, g^,^{x,u) is algebraically determined in terms of g/^^, and the second 
term in (1) describes just a potential for the induced metric on the brane. 

Of course the choice of this second boundary condition is by no means unique 
at the level of auxiliary extra dimension: while g^i,{x,u = 1) = 7]^^ was originally 
chosen to describe a particular completion of Fierz-Pauli massive gravity, nothing 
prevents us from considering a more general boundary condition to accommodate a 
larger family of completions. In this way we have written the most generic gravi- 
ton's potential in terms of a geometrical construct that naturally arises in higher 
dimensional theories of gravity, where /c^j, is taken to be the extrinsic curvature. 

The abovementioned generalization of the boundary condition is straightforward 
and will be the subject of §2, where we show that one can consistently take g^ui^, u — 
1) a generic function of g^^, the 4D-brane metric. In §3 we perform the dimensional 
reduction to obtain the 4D potential for graviton, and explicitly show that in the 
original choice of boundary condition the absence of the ghost does not persist at 
higher than cubic order (in agreement with [14]). Thereafter a recipe will be provided 
to reproduce the potential of [3] by a proper choice of the boundary condition, 
ensuring the stability of the full theory in the decoupling limit. In §4 we consider 
the generalization of the model to higher extra dimensions and prove that in the 
most tractable case with a spherically symmetric boundary condition in the bulk. 
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one and many extra dimensions are equivalent. 

We further provide a simple explanation for the absence of the cubic ghost 
and the presence of higher order ones in the original choice of boundary condi- 
tion {g^i,{x, u — 1) — r)ni,), in appendix A. In appendix B we address the possibility 
of modifying the bulk action by inclusion of other geometrical constructs so as to 
avoid the adjustment of the boundary condition. It is shown there that the most 
natural candidate which is (the auxiliary version of) Gauss-Bonnet term, as the only 
"Lovelock" term other than Einstein-Hilbert in 5D for which the variation principle 
is well defined, cannot do much better than the original action in predicting right 
coefficients. Nevertheless one can build a ghost-free theory by adding higher powers 
of kfj^i, with specially tuned coefficients. Finally, in appendix C we give the most 
general boundary condition which gives healthy theory at quartic order. 



2 Generalized Boundary Conditions 

We begin this letter by considering the generalization of the boundary conditions for 
(1). For convenience we will assume the location of the boundaries to be at w = 0, 1 
rather than u = ±1. As usual we fix g^ij{x^u = 0) = g^v{x), while the choice 
at M = 1 is in principle arbitrary since this boundary does not have an intrinsic 
dynamics and the values of the fields on it are completely determined by imposed 
conditions. Hence, the most general boundary conditions would be 



g^,y{x, M = 0) = g^y{x) and g^y{x, u = l)= g^u{gaf}), 



(2) 



where g^j^iy — 7]^^ has to vanish when g^^ = 1]^,,, if the theory is to describe massive 
gravity around Minkowski. The variation of the action (1) then generalizes to 
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gG,Jgnu=o-m' ^-~g{k>^'^ -~g^-k)5g 



ip,v 



u=0 



du~sj —g 



^vhfxv g^i/duk 'T^Qiivik^^ -\- k ) 2k'^kai> + kkfj,i^ 



S~g 



from which it follows that the bulk [u > 0) equation of motion is 



r9 k 



kkn 



the same as in the case with field independent boundary condition at u 
brane equation of motion on the other hand gets slightly modified to 



IJ,V 



m 



dgai3 
1 dgp^ 



(3) 

(4) 
1. The 

= 0. 



Notice that this differs from the 4D equation of [10] by the last term, which arises 
because of the non- vanishing Sgfj,i,\u=i^ ■ 



^It should be mentioned that this does not make the variation principle ill-defined. 
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An equivalent and often simpler way of analyzing the AED models is to integrate 
over the w-coordinate using the solution to the bulk equation (4) with the boundary 
condition (2) in order to obtain an effective 4-dimensional potential for graviton. 
This was the strategy pursued in [13] for the special case of g^vix, u — 1) — r]^i,, 
treating h^i, perturbatively, but the generalization is straightforward. One simply 
needs to write 

9niy = Vixu + = g^^ + h^^ , with /i^^ = r^^^ + /i^^ - g^^, (5) 

and expand all g^^, around g^,^ so that all computations go through identically if the 
following replacements are made: 

Vfiu Qtiv , and h^,^ h^^ . (6) 

Therefore one just needs to take the original 4D effective potential and make the 
substitution (6) (including the overall factor sj— det(77) = 1 — >■ •>/— det(^)) . 

3 Quartic Ghost and a Remedy 

Perhaps the simplest method to check models of massive gravity against Boulware- 
Deser ghost [1] is to take the routes of refs. [13] and [15], namely, to write the 
graviton potential in the unitary gauge, expand it around the Minkowski background 
and compare it with the unique set of conceivably ghost-free potentials constructed 
in [3]. Obviously to perform this comparison the latter should also be re-expanded 
in the unitary gauge where no Stiickelberg fields are present. 

Following the logic of the previous section, in order to perform the dimensional 
reduction perturbatively in h^y, we define 

li) = g^y{x, u) - g^y = Hj^J{x, u) + Hj^J(x, u) + . . . , (7) 

where Hjjl^ is of n*^ order in h^u and except Hj^J which is given by 

Hl^J{x,u)^{l-u)h,y{x), (8) 

all the rest vanish both at u = and m = 1, so that g^u{x,u) satisfies the desired 
boundary conditions at each order. This property of higher order bulk solutions 
together with the absence of second or higher derivatives of fields in the bulk action 

ensures that to obtain the n**^ order graviton potential one only needs to solve h^^, 
up to (n-2)'^'^order. To wit note that the action (1) starts from quadratic in h^i, 
which thus contains the n*^'^ order term of the schematic form 

du duHj^^-'^d^Hj^J = Hj--'^duH;^J\l - f du H;--''>dlHj^J , (9) 
Jo Jo 
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but the first term on the r.h.s. vanishes because of the vanishing of on 
the boundaries, and the second term because of the linear equation of motion which 
is satisfied by H^p^J{u). 

In the case of field-independent boundary value g^i, — 7)^^, the comparison with 
the tuned polynomial was carried out in [13] up to cubic order yielding agreement 
with one class of (decoupling limit) ghost-free potentials characterized by C3 = 1/4. 
Here we show that this agreement breaks down at quartic order and in the appendix 
A using a simple covariantization we explain why it could have been expected. We 
also find the most general modification of boundary condition that give rise to the 
ghost free 4D action up to quartic order. This is achieved by perturbative adjustment 
of the boundary value at -u = 1, hence one can continue this process to all orders. 

As already mentioned, to obtain fourth order 4D effective potential one needs 
only up to the second order bulk solution: 



12 



VnAKp -h)- h^K^ + -hh 



(10) 



which after some work yields the following unitary-gauge potential 




.Ml - VTjChLf + ^h'hl - :^h') , (11) 
12 144 72 144 / 

with indices contracted by matrix g^^", inverse to Qui,. This expression should be 
compared to 

VtuLiH,.) = V^[Hl-H' + c,Hl + C2HHl + csH' 

+d,Hi^ + d,HHl + d,{Hl^f + d.H^Hl^ + d,H^] , (12) 

where H^^, which is a covariant version of h^i, = g^u — rj/^,^, reduces to h^i, in 
the unitary gauge, indices are now contracted with the full metric g'^''' and all the 
coefficients are expressed in terms of C3 and d^ as given in [3]. 

Concentrating on the case of constant boundary g^j^^ = 7]^^^ where we know that 
C3 = 1/4 guaranties up to 3''^order agreement between (11) and (12), and, expanding 
the latter in the unitary gauge around Minkowski background we get 



tuned I 'Vf' J — 'Vf ''' '^ixv^ ^'^'^ixv ^'^ 



+dihl^ + d2hhl^ + daihl^f + d^h^hl^ + d^U^ . (13) 

However this can never coincide with (11) for any value of 0^5 and consequently 
the theory suffers from a ghost at quartic order. This can most easily be seen 
by introducing Stiickelberg fields in terms of which there remains a quartic self- 
interaction of the helicity-0 part tt of the schematic form 

Ml,m\ddi.t = j^.iddn^r ^ ^i^dn^f ' (14) 
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where tTc = Mpim^Tr = AgTr is the canonically normahzed field. In terms of interac- 



ts' 



tions of TT the agreement with the adjusted potential to 3 order translates into the 
absence of the cubic interactions: 

-{ddn^f ^ ^Jddn^r , (15) 



which are generically present in theories of massive gravity [16]. Consequently (14) 
is the most strongly coupled interaction in the theory or equivalently the scale 
A4 = (Mpim^)^/^ is the smallest mass scale by which any interaction in this model 
may be suppressed. We can therefore define a decoupling limit 

Mpi ^- 00, m^O, A4 - fixed, (16) 

in which (14) is the only interaction that survives and it propagates ghosts lighter 
than the cutoff around any reasonable astrophysical background (as shown in [2] for 
the similar case of the cubic interaction (15)). The fact that the action contains 
only a finite number of terms leaves no ambiguity in the instability of the model 
since no resummation may be invoked to remove the ghost. 

Having discussed the problems of the original model let us comment on (11) with 
general boundary condition. A similar analysis shows that ghost can be avoided by 
order-by-order adjustment of the boundary condition 



+ hfh,X -rhf[h\.h-rhf\h\. 



(17) 



The requirement of the Fierz-Pauh structure constrains hi^2-i the absence of the cubic 
Boulware-Deser ghost relates the quadratic coefficients &i^.^..^4 to each other and so 
on. The most general healthy coefficients of the above expansion are explicitly given 
in appendix C. Here, on the other hand, we give just one of the simplest expressions 
for gf^i/ which cures the instability of the theory in quartic order 

%v = r]^,^ + -^V^.uhlp - ^Vi^uhhlf^ + ^r]^^h? - ^h^^hl^ + (18) 

The theory with this choice of the boundary condition, once reduced to 4D gives 
the 4*^^ order potential of [3] corresponding to C3 = 1/4 and d.s = 0. 

We would like to stress that our approach is perturbative in contrast to [14], 
where authors have performed the w-integral exactly and obtained nonlinear 4D 
action as a function of Q/^i, (in their notation the metric at = 1 is labeled as 
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f^u)- The advantage of that framework is that one may try to find an exact ghost- 
free boundary condition by equating F{g'^'^g^a) of [14] to the ghost-free potential 
^gh-iridtMu, Vfiu) of [4]. However, because of the transcendental nature of the equation 

[det(^?^"^„.)]^/' - 2 [detig^'^ga.)^' x 

cosh ^jTiMg^g^,)]^ - ^[Tr \n{g^gau)?^ + 1 = ^Wgh-fr, (19) 

one seems to be forced to do the perturbative analysis. 



4 Generalization to Multi-D 

Another natural generalization of AED is to consider a multi-dimensional auxiliary 
space instead of a single dimensional one. In this section we investigate the special 
case when there is a rotationally invariant condition on the boundaries of the d- 
dimensional bulk, and show that it gives identical 4D effective action as the original 
AED model. 

Using the cartesian coordinates ya, where a = 1,2, ... ,d, the action of the extra 
dimensions generalizes to 

S,= I d'y^g{kl^-kl), (20) 

where kafj,u = \dag^Mu and summation on repeated indices is implied. We are willing to 
impose a spherically symmetric boundary condition, however, one should be cautious 
in this case since the solutions of the Laplace equation are singular at the origin. 
This is a generic feature of frameworks with co-dimension >1 branes [?], and it is 
well known that the singularity can be regularized by assigning a finite width e to 
the brane. It suffices to take the radial coordinate, defined as r = ^JyaVa-, to range 
in the interval r G [e, 1]. Therefore the region r < e is excluded from the integral in 
(20), and the boundary conditions are modified to 

~9p.v{x, r = e) = g^^, g^vix, r^l)^g^y. (21) 

From this boundary condition it follows that the solution to the bulk equations 
of motion will also be spherically symmetric and therefore (20) can be simplified 
considerably by going to the spherical coordinates. Integrating over angular variables 
and dropping an unimportant normalization constant which can always be absorbed 
in the definition of the graviton's mass we obtain 

Sd = J' r^-'dr^g [{K,,f - kl] , (22) 
with kruv = \drgiji,v However by a change of the variable of integration to 

d-=^.^ (23) 
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and a rescaling such that u ranges in [0,1], the action (22) transforms to that of 
single extra dimension model (1). Hence the spherically symmetric multi-D analog is 
equivalent to the original model with one extra dimension. However, this equivalence 
does not generically persist the modifications of the bulk action, see appendix B for 
more details. 
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Appendix A. Covariantization by A/^ 

As first pointed out in [11], AED models with fixed (e.g. rj^i,) boundary condition 
at M = 1 have their own natural candidate for playing the role of Stiickelberg fields, 
namely the ADM shift vector iV^. Consider Einstein- Hilbert action in 5 dimension 
with metric Gmn and decompose it on spatial slices of constant 5**^ coordinate using 
ADM parameters: lapse N = (G^^)"^/^, shift Nf^ = Gs^, and 4D metric Qf^^ = G^i/ 

S = 2M2 J d'xdu^-det{g^,)N[R^^^ {g) - K^^'K^, + K^] , (A-I) 



where 



K^. = 2^(5«^M- - D^N, - D,N^) , (A-H) 



all indices are raised using {g'^Y" = g'^'^ , and is the covariant derivative with 
respect to the 4D metric. The action is invariant under the full 5D diffeomorphisms 
and in particular the subclass of li-dependent 4D diffeomorphisms 

u^u' = u, ^ x"'{x, u) , (A-111) 

under which N and R^^^ behave as scalars which implies that i^T^jy must transform 
as a covariant tensor. 

Likewise, we can restore u-dependent diffeomorphism by replacing k^, in the 
action of AED with Kn^, — (9u^^^/2 — D(^^Ni,)) and stipulating that N^, which is 
now some auxiliary field, transforms the same way as the 5D gravity's shift vector 
does, namely 

N" —— + d^x'" . ( A-IV) 

ox'^ 

Having restored this class of diffeomorphisms the action is now invariant if one 
reparametrizes the 4D metric g^iy on ti = brane but keep the other boundary fixed 
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at g^u{x,u = 1) = rj^y. In other words covariantizes the 4D effective Lagrangian, 
and can be regarded as the Stiickelberg field. This covariantization fails to work 
in the more general case where the boundary condition aX, u — 1 {g^iy) depends on 
g^ti, because now the reparametrization at w = brane changes the g^^, but not 
necessarily in the manner of a coordinate transformation. 

In terms of it is easy to understand the absence of the cubic ghost and the 
emergence of the quartic one: Working around Minkowski where A'^^ is also taken 
to be small, the first order bulk equations of motion become 

5X.-29„9(^iV,) = 0, (A-V) 
dadu(hl - S;h) - d'^id^N^ - d^N,) = , (A- VI) 

which are solved by h^jjil = (1 — u)h^j,i, (as in the non-covariant case (8)), and a 
constant along the u direction. This linear solution, as before, is sufficient to 
obtain cubic effective action which consequently contains at most two powers of 
N^, because there are originally two A^^u's present in the action and lij^} does not 
contain any. This explains the absence of the cubic ghost since there cannot be any 
cubic self-interaction of the helicity-0 mode which is contained completely in the 
Stiickelberg field A"^ in the decoupling limit [11]. This argument, however, breaks 
down beyond that order because higher order solutions will necessarily contain A^^j 
(e.g. }ij}u D A"^^) which upon reduction to 4D result in quartic and higher order 
terms in A"^, explaining the presence of the interaction (14). 

Appendix B. AED with 5D Gauss-Bonnet Term 

In this section we extend the framework of the auxiliary extra dimension (AED) by 
addition of terms descendent from the Gauss-Bonnet (GB) action. The latter, being 
the only other five-dimensional Lovelock invariant besides the Einstein-Hilbert term, 
is given by 

Cqb — {R^ — '^RabR^^ + RabcdR^^^^) , (B-I) 

with yl, . . . = 0, 1, 2, 3, 5 and n being an arbitrary constant. 

In order to find the AED analog of (B-I) we impose 51^5 = and g^^ = 1 on the 
metric tensor. Furthermore, since we are not interested in generating derivative self- 
interactions for graviton, we remove terms containing four-dimensional derivatives. 
As a result (B-I) reduces to 

^gT = ^g^'d.k,, {k' - kip) + 2d,k,, (n'^ - k'^^k) 

+ \ {-Ukl + l&kkl, + IklXp - lOkXu + k')] . (B-II) 

However, from the definition of A:^,^ it follows that there are terms with more than 
one derivative per field in (B-II), hence, one needs to introduce a boundary term 
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in order for the variation principle to be well-defined. Including those we get the 
following modification to the Einstein-Hilbert Lagrangian 



K 



3 



n=l 



T7^(2<. - 3A;A;^, + k') - / du^-~g{kl - k' + C 



u=0 Jo 



^/ii/ ^ ^GB 



(B-III) 



In order to integrate out dimension, one has to find the solution to the bulk 
equations of motion which now generalizes to 

QiJ.i'^uk — duk^n — ~-^Qfj.v{k + kg^^) — ^kfj^k^^ + kkfj^u 

\^yk^i^k ~\~ d^ki^^k ki^^ ~\~ k^^^k k^^^k j ^ , (B-IV) 

where we have presented the contribution of the GB term schematically since it is 
sufficient for our purposes, as will be seen shortly. We choose the boundary condition 
to be g^^uix, u^l) = r]^^. 

The easiest way of solving (B-IV) for h^„{x,u) = gf^,y{x,u) — rj^i, is to proceed 
order-by-order in four-dimensional metric perturbations h^y{x) = h^y{x^u = 0). 
One immediately notices that the newly added terms proportional to k in (B-IV) 
start to change the solution only at the 4*'* order in h^y{x). Using therefore the old 
linear solution (8) it follows that the Gauss-Bonnet term does not contribute to the 
cubic 4D effective action: The cubic bulk terms in (B-II) lead to dl^H^^^ = 0, while 
the boundary terms in (B-III) evaluated to S'^'^order are identical at ti = and u = 1 
and cancel each other. 

As in §3, to find the 4D action up to 4*'* order one only needs the 2"^order bulk 
solution (10) which leads to the following four-dimensional Lagrangian 

r2 r^U "^^^ll r^.. fu2 u2 u3 , ^UU2 ^ u3 



L = Mpiv^it:-^^v^x(/i^,-/i^-/i^, + -M^,--/i^ 

+^ (6(22 + 3«;)/iJ, - (53 + ^n)hl,hl^ - 12(11 + 2K)hh% 

+2(29 + 9K)h^hl^ - (5 + 3K)h^) + 0{h^)), (B-V) 

with indices contracted by the Minkowsi metric rj^'^. It is easy to see that this will 
never match the tuned potential of [3] for any value of k, meaning that the quartic 
ghost-like pathology of the original model can not be cured by terms of geometrical 
origin, since the GB term is the only available one in 5D. 

So far we limited ourselves to potentials which are motivated by some 5D geomet- 
rical construct, however by giving up that requirement one can naturally generalize 
the potential term in (1) to 

r+l 

V{g^,.) = mM duy^ {kl, - A;' + ai/cj, + aa/c/cj, + a^k^ + ...), (B-VI) 
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where now the coefficients at each order should be chosen such that after reduction 
to 4D and introduction of Stiickelberg fields pure helicity-0 interactions add up to 
a total derivative at each order. 

One interesting observation that follows from the covariantization (appendix 
A) is that the n**^ order terms in k^y make no contribution to the (n+1)*'^ order 
ghost-like interactions. This is because as was the case for the cubic effective action 
in the original model, only the first order solution is needed to be substituted in 
{kiiu)^ terms. After covariantization (i.e. replacing k^i, with K^i,), HjlJ remains 
independent of N/^ and therefore the highest power of N^i at (n+l)^*order is (A^/^)". 
This, however, cannot possibly affect (n+l)***order self-interaction of helicity-0 mode. 

It is also worth mentioning that the equivalence between one and several spher- 
ically symmetric auxiliary dimensions does not survive general modifications of the 
bulk action that include higher powers of k/j^i,. The addition of k^^^, terms causes 
the multi-dimensional model to deviate from its co-dimension one counterpart at 
(n+l)^*order in perturbations. Similarly the Gauss-Bonnet terms lift this degener- 
acy because of containing higher powers of dr, moreover there is one new GB term 
for each extra dimension which can in principle be included in the action. 



Appendix C. Fine- Tuning of the Boundary 

There are two sets of coefficients (17) that do not give rise to the ghost. One of 
them is given by 

b? = 0, v4Vi, 

b? = = f + l(l-4c3-26«+4c36« + 6«^), V^r, 

6f = l(6«-l)(-l + 4c3 + 26«). 



b? = ^ + ^(8c^(6(^'-l)-(6«-l)(-3 + 164 + 36«)-4(3 + 46«)6f 



+4c3(4+(-5 + 6^^0^2 +46r)) 

b?^ = —± - + — ((6^^^-l)(-96cl + 4324 + 12c3(27-56«) 

18(6^^ - 1) 864 

+{b^^^ - 1)(61 + 56^'^)) + 24(9 - 10c3 + 166^'^)6f + 288(6^ - fef )) 

^(3) ;f2)2 

+12c3(3 + 6«) + (ftW _ i)(7 + 56^)) + 48(C3 + 2&«)6f ) 
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&f = , ' ,(1), + ^((b? - imcl + 1084 + 12c^b? + {b? - 1)(1 + 2b?)) 
3{-l + b?) 72 

+66f (-9 + 2003 + 46^)), 
4'^ = ^^^(7-724-146(^^+4(-3c3(3 + C3) + 6c36^') + 6(^)')), V 6g. 
while the other one being 

b?' - -f -%^(-8 + 16c3 + 36«), V6r\ 
fe?^ = ^((-l + 2c3)(-l + 6^^^)-26f), 

&f = ^(-l + 4c3 + 26(^)), 

6f) = :^(4(6^'' - l){cl + 224) - 4c3(13 + b?i-16 + 3b?) + ib?) 
+6^'V23 + b?{5 + Ab?) + m?) + 2(7 + 66^ - 86^)), 



1,(2) 



2 



^'f = --^-^ + -^((l-6?))(-85 + 408c^ + 10804 + 12c3(27 + 76«) 

18(-l + 6^0 ^^"^ 
+6^') (14 + 536^'^)) - 12(-9 + 40c3 + 20b?)b? - 144(26f + ftf )), 

7(2)2 

6f = , ' ,(1), + ^((4^^ - l)(48c^ + 1444 + 12C3(3 + b?) 
18(-l + 6^^) 8o4 

+{b? - 1)(7 + 5b?)) + 24b?{2c, + b?) - 1446f ), 

2l{2)2 

6f = — ^-^ + -((4')-l)(84ci + 1084 + 12c36« + (6«-l)(l + 26(^))) 
3(_l + 4i)) 72 

+66^ (-9 + 20c3 + 46^'))), 
b? = ^^2r^(7- 724- 146^') +4(-3c3(3 + C3) + 6c36^^) + V fcg. 

These coefficients give the most general boundary condition at u = 1, for which the 
theory is ghost-free (in decoupling hmit) up-to 4*^ order. It is quite straightforward 
to continue this tuning to arbitrary order. 
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